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We apply the coordinate-space method by Liischer and Weisz to the computation of two-loop diagrams in full 
QCD with Wilson fermions on the lattice. The essential ingredient is the high-precision determination of mixed 
fermionic-bosonic propagators. 



At the Lattice conference of last year we pre- 
sented an algebraic algorithm that allows to 
express every one-loop lattice integral with gluon 
or Wilson-fermion propagator in terms of a small 
number of basic constants which can be computed 
with arbitrary high precision Q . This was a gen- 
eralization of what we did previously with purely 
bosonic integrals || and it was also an essen- 
tial step in order to apply the coordinate-space 
method by Liischer and Weisz Q to higher-loop 
integrals with fermions. 

Let us consider, in order to fix the ideas, an 
example of a two-loop integral at zero external 
momentum, like 



where 
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trapolate to infinite L. A possibility is rewriting 
@) as 

1 = T8 /, J^o D F {l)D F {r)D F {l + r) (3) 

where each component 1^ and Ti runs over the set 
27r(n + 1/2) /Z, n = 0, 1. From the sum 

we exclude the points such that I + r = mod 
2-7T where the third propagator diverges. Another 
possibility is to use (j^) but with li and r.; running 
over the set 2im/L, n = 0, . . . , L — 1. In this case 
however there are more problems with the zero 
modes and one should exclude from the sum the 
terms with I — r = I + r = mod 2ir. For this 
reason we have decided to use the first method. 
For our previous example we get for increasing 
values of L 



defines the usual Wilson propagator with momen- 
tum I at zero fermion mass with Wilson parame- 
ter rw ■ In all our calculation we have set rw = 1 
but the method can be applied for any value of 
rw- A possible strategy to evaluate such an in- 
tegral amounts to replace each integration with a 
discrete sum over L points and afterwards to ex- 



L = 10 
L = 18 
L = 20 
L = 26 



0.000799652 
0.000848862 
0.000853822 
0.000863064 



(4) 



Then, using an extrapolation of the form 

ai log L + a 2 a 3 log L + a 4 
«o + T2 + Z 3 

and data with 6 < L < 26, we obtain the estimate 
/ w 0.000879776 



2 



which has to be compared with what we obtained 
by using the coordinate-space method 



0.0008797779181(12) 
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where £>s(fc) = k is the usual bosonic propa- 
gator on the lattice. In the coordinate-space ap- 
proach we are interested in the evaluation of the 
lattice sums 

I(p 1 ,q 1 ,a,p 2 ,q 2 ,b,p 3 ,q 3 ,c) = (6) 
G{pi,qi,x + a)G(p 2l q2,x + b)G(p 3l q 3 , X + S) 

X 

In this notation our previous example corre- 
sponds to 



7 = 1(1,0,0,1,0,0,1,0,0) 



(7) 



In the evaluations of these sums we make use of 
the following advantages: 

• only four infinite lattice sums must be com- 
puted; 

• the G(p, q, x)'s can be determined with the 
desired precision, for a sufficiently large do- 
main of values of x, by using our algebraic 
algorithm 

• the asymptotic expansion for large values 
of |x| of the G(p, q, x)'s is easily computed. 
For example 



G(l,0,a?) = 4 



1 

4x 2 



1 



2X4 



10 



52 

"^5 



160x1 



192x 6 



(8) 



where x n = 



Let us now consider how to compute sums of the 
type 



(9) 



on the lattice A. Of course we will not be able 
to sum over all the lattice. If |x|i = ^ |a:^|, 
we will perform a sum over a domain of the type 
D p = {x G A : |x|i < p}. The problem is to give 
an estimate of the error. If f(x) decreases for 
large \x\ as l/|x| 2fc we expect the sum restricted 
to D p to behave as 
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Thus we will estimate 
P 



2k -A 



S(p)-E(p-1) 



(10) 



(11) 



Our error formula seems to work correctly. We 
will use this formula to estimate the error on the 
integration sums. We can also define an improved 
estimate for S by 



S(p) 



P 



2k 



r (E(p)-E(p-l)) 



(12) 



Notice that now our error estimate is very conser- 
vative. The larger is k the best is the estimate. 
For this reason, if we know the asymptotic be- 
haviour Af(x) of the function /(x) for large x it 
is convenient to write 



E«£[/(x)-4f(z)]+£Af(*) 



(13) 



because the difference is decreasing faster at in- 
finity and the sum of Af(x) can be computed 
directly on the infinite lattice by using harmonic 
polynomials and ^-functions as explained in Q. 
Coming back to our preferred example, by sub- 
tracting an increasing number of terms of the 
asymptotic expansion, we get the estimates 

I = 0.0008798104043 ±0.0000008730034 

= 0.0008797776858 ± 0.0000000029778 

= 0.0008797779227 ± 0.0000000000410 

= 0.0008797779181 ± 0.0000000000012 

In our work we have always used the asymptotic 
expansions to increase the precision of the es- 
timates. In each case we have subtracted the 
asymptotic behaviour to order l/|x| 10 : therefore 
the function which is summed over a finite lattice 
decays at least as l/|x| 12 . 



3 



A number of checks have been performed on the 
table of numerical integrals that we have collected 
(at the moment we have a number of entries of 
order 10 4 ). In particular 

• In the case in which we restrict ourselves 
to purely bosonic integrals we compare per- 
fectly with the numbers given in B . 

• Because of translation invariance, for every 
v on the lattice 

I(p 1 ,q 1 ,a,p 2 ,q 2 ,b,p 3 ,q 3 ,c) = (14) 
I(pi,qi ,a + v,p 2 ,q2,b + v,p 3 ,q3,c + v) 

• From the definition of the bosonic propaga- 
tor one easily gets 
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which can be used to derive relations among 
different Fs. Similarly, other relations 
can be obtained from the definition of the 
fcrmionic propagator. 

• New identities are obtained by integration 
by parts, that is by using 



(16) 



All the checks we have performed are satisfied 
with a precision of at least 10 -10 

We have prepared a completely automatic pro- 
cedure which evaluates Feynman diagrams at two 
loops. It goes through the following steps 

1. Each diagram is reduced as a sum of the 
integrals /'s, previously defined in (^J). A 
mass is added at all propagators in order to 
regularize the infrared divergences. 

2. All the possible symmetries (cubic: transla- 
tions, permutations of the axes, inversions 
of the axes; permutations of the three prop- 
agators) are used to reduce the number of 
terms. 



3. Through subtractions all the terms are writ- 
ten as convergent sums plus product of 1- 
loop integrals. 

4. All 1-loop integrals are expressed as in 

5. The convergent sums are replaced by their 
numerical estimate obtained from a pre- 
compiled table of lattice sums in the domain 

^21- 

This procedure is now being used to compute the 
mixing coefficients of the four fermion operators 
of the lattice weak hamiltonian with the dimen- 
sion 5 operators. We have chosen this compu- 
tation because it starts at the two-loop level and 
because these coefficients have already been stud- 
ied by using the momentum-space approach 
The computation is highly non trivial. The typ- 
ical input, for each diagram, contains order 10 4 
terms and produces a final result of order 1CT 4 . 
Unfortunately at the moment our results are still 
preliminary and we are performing all possible 
checks on our evaluations. 
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